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ABSTRACT 



The way galaxies assemble their mass to form the well-defined Hubble sequence is amongst the most debated topic in modern 
cosmology. One difficulty is to link distant galaxies, which emitted their light several Gyrs ago, to those at present epoch. Such a link 
is affected by the evolution or the transformation of galaxies, as well as by numerous selection and observational biases. We aim to 
describe the galaxies of the Hubble sequence, 6 Gyrs ago. We intend to derive a past Hubble sequence that can be causally linked to 
the present-day one, and further estimate the uncertainties of that method. We selected samples of nearby galaxies from the SDSS and 
of distant galaxies from the GOODS survey. We verified that each sample is representative of galaxies selected by a single criterion, 
e.g., Mj(AB) < -20.3. We further showed that the observational conditions needed to retrieve their morphological classification are 
similar in an unbiased way. Morphological analysis was done in an identical way for all galaxies in the two samples. We found that 
our single criterion is particularly appropriate to relating distant and nearby galaxies, either if gas is transformed to stars in relatively 
isolated galaxies or, alternatively, if they accrete significant amounts of gas from the intergalactic medium. Subsequent mergers during 
the elapsed 6 Gyrs, as well as evolution of the stellar populations, are found to marginally affect the link between the past and the 
present Hubble sequence. Indeed, uncertainties from the above are below the errors due to the Poisson number statistics. We do find 
an absence of number evolution for elliptical and lenticular galaxies, which strikingly contrasts with the strong evolution of spiral and 
peculiar galaxies. Spiral galaxies were 2.3 times less abundant in the past, which is compensated exactly by the strong decrease by 
a factor 5 of peculiar galaxies. It strongly suggests that more than half of the present-day spirals had peculiar morphologies, 6 Gyrs 
ago, and this has to be taken into account by any scenario of galactic disk evolution and formation. The past Hubble sequence can be 
used to test these scenarios and to test evolution of fundamental planes for spirals and bulges. 

Key words. Galaxies: evolution - Galaxy: formation - Galaxies: spiral - Galaxies: peculiar - Galaxies: starburst - Galaxies: bulges 



1. Introduction 

Galaxies are complex objects containing several tens of bil- 
lions stars, as well as gas and dust. Present-day massive galaxies 
are regular and relaxed systems, and are made of a dispersion- 
supported bulge that is surrounded by a rotatio nally supporte d 
disk. They all fit the initial scheme proposed by [Hubble! (Il926l) . 
which is called the Hubble sequence. In the Hubble scheme, 
which is based on the optical appearance of galaxy images 
on photographic plates, galaxies are divided into three general 
classes: ellipticals (E), lenticulars (SO), and spirals. In the se- 
quence, galaxies are organised from pure bulge (the elliptical) 
to increasing disk contribution to their light or mass (SO to Sc). 
Further improvement has defined a fork that splits barred (SB) 
from non-barred spirals. There are also few massive galaxies, 
less than 10%, showing irregular or peculiar morphologies, and 
these escape to this classification. 

How was the Hubble sequence in the past? The Hubble 
Space Telescope (HST) revealed the morphologies of distant 
galaxies, and a simple inspection of the deepest observed fields 
show s how galaxy morphologies have evolved (van den Bergh, 
l200l . Indeed, we do not lack morphological details as the 
HST/ ACS instrument allows investigation of the details of 200 
pc at z ~ 1 ■ However, this does not suffice to establish the past 
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Hubble sequence. First, the cosmological dimming is so severe 
that only long exposures, such as in GOODS or UDF, can de- 
tect the optical radius of a Milky Way like galaxy at z=0.4-0.5 
and 1.3-1.5, respectively. Second, the past Hubble sequence has 
to be made of galaxies that can be connected with the present- 
day ones: sporadic events such as galaxy collisions could have 
dramatically affected this link. Third, although distant galaxies 
are linked to present-day ones through the isotropic principle, it 
needs a simple (and preferentially single) criterion in selecting 
them in two different redshift slices. 



The motivation of this study is to establish the Hubble se- 
quence in the past and to verify to what accuracy it can be 
linked to the sequense at the present epoch. A considerable ad- 
vantage would be to compare local galaxies to those at moderate 
redshifts, in order to limit most of the complexities discussed 
above. At z=0.65, corresponding to an epoch of 6 Gyrs ago, the 
GOODS imagery -after correction for cosmological dimming- is 
as deep as or deeper than the SDSS imagery at z=0. Moreover, 
it is unlikely that major mergers could significantly affect the 
number density of galaxies since that epoch. The selection cri- 
terion has to be kept simple and should address the mass in one 
way or another. There is a relatively good correspondence be- 
tween stellar mass and absolute magnitude in the near-IR, so we 
choose a limit Mj(AB) < -20.3. Such a value is adopted from 
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the IMAGESQ survey and corresponds approximately to M ste u ar 
> 1.5 10 10 M o . However, one needs to make it clear whether this 
simple criterion select similar galaxy masses in the two redshift 
bins. 

The methodology for morphological classification is a cru- 
cial issue, especiall y when comparing galaxies at two differ- 
ent epochs. Quoting [Kormendv & Kennicuttl (120041) . "it is use- 
ful only if classification bins at least partly isolate unique 
physics or order galaxies by physically relevant parameters. 
The Hubble-Sandage-de Vaucouleurs classification scheme has 
done these thing s remarkably well". We thus adopt a method- 
ology (see, e.g., iNeichel et al.L 12.0081) that systematically com- 
pares morphologies of distant galaxies to those of local galaxies. 
Identification of peculiarities would thus be directly assessed by 
the discrepancy between a given galaxy and the galaxies that 
populate the local Hubble sequence. Such a technique is based 
on a simple and reproducible decision tree and a full decompo- 
sition of the galaxy light profile in tw o components , name ly the 
bulge and the disk. A striking result of Neic hel et al.l d2008) is the 
excellent agreement between such a morphological classification 
and classification of their kinematics. By studying a representa- 
tive sample of 52 distant galaxies, they demonstrate that 80% of 
rotating galaxies possess spiral-like morphologies, and that 95% 
of galaxies with complex or chaotic kinematics have peculiar 
morphologies. In other words, pe culiar morphologies are asso- 
ciated with anomalous kinematics. [Neichel et al.ld2008l) has also 
shown that other classifications based on automatic assessment 
of compactness or asymmetry do not correctly predict the kine- 
matics, as they overestimate the number of spiral- like morpholo- 

[ ies. Such an effect h as indeed been notic ed by IConselice et al.l 
20031 120051) (see also lKassin etafll2007l) . 

Finally, when comparing galaxy morphology at different red- 
shifts, one must pay attention to all possible biases linked to the 
use of different instruments and, more importantly, must min- 
imise morphological ^-corrections. In the present paper, after 
considering the above effec ts, we general i zed th e morphological 
classification method from Neich elet al.1 ([2008) and we applied 
it to a local sample of galaxies observed by SDSS and to a distant 
sample of galaxies observed by HST/ACS, in order to compare 
their properties. 

Samples and their selection are described in Sect. 2. 
Methodology and data analyses are shown in Sect. 3. Results 
of our morphological classification are presented in Sect. 4, and 
Sect. 5 discusses to which accuracy one may establish the past 
Hubble sequence. Finally, in Sect. 6 we summarize our conclu- 
sions. Throughout this paper we adopt cosmological parameters 
with Hq = lOkms^Mpc' 1 , Qa = 0.7. Magnitudes are given in 
the AB system. 



high resolution images in at least three optical bands. These fur- 
ther limitations are only rejecting objects for technical reason^ 
in the local Universe since data from SDSS have reached an in- 
comparable homogeneity. While the GOODS image quality and 
depth are as well homogenous for distant galaxies, spectroscopic 
surveys s how quite significant h eterogeneities. We followed the 
steps of Ravikumar et al.l (120071) showing the ability in gather- 
ing a representative sample from the collection of spectroscopic 
samples in the GOODS area. All in all, the luminosity functions 
of both samples are compared to those of similar surveys, and an 
excellent agreement is found. 



2.1. Local galaxy sample 



Galaxies were selected from the Fukugita et al] d2007l) sam- 
ple. It includes 2253 galaxies with Petrosian magnitude lower 
than rp=16, that are lying in a 230 deg 2 rectangular area in the 
equatorial area of the northern sky. From this sample we only 
kept galaxies with availabl e spectra and imag ing from the Sloan 
Digital Sky Survey (SDSS:lYork et al.U 2000) Data Release three 
(DR3; Abazajianet_aL|, I2005I) . All the objects have images in 
u (3551 A), g (4686 A), r (6165 A), i (7481 A) and z (8931 
A) bands. Relative ma gnitudes were retrieved from the 2MASS 
dSkrutskie et al.L 120061) and SDSS-DR3 catalogs. It results a se- 
lection of 2113 galaxies. We then considered only the 1665 
galaxies with an absolute magnitude My < -20.3, following our 
selection criterion. Figure 1 shows the My ys redsh ift plane of 
the local SDSS catalog from Fukugit a et al.l d2007h . There is a 
clear lack of intrinsically faint galaxies at z > 0.03, especially 
blue galaxies with magnitude just below rp=16. This is illus- 
trated by the blue curve which corresponds to a flat spectrum 
between r-band and J-band, i.e., that of a pure starburst at the 
sample magnitude limit. Figure 1 shows that at z < 0.03, all 
My < -20.3 are included in the Fukugita et al. (2007) sample. 
It represents a sample of 218 galaxies (red points). Finally, from 
these 218 galaxies, we only kept those whose spectra include 
the [On]/13727 emission line, further reducing the sample to 1 16 
galaxies. These galaxies have a redshift of 0.0207< z <0.030. 

Figure 2(a) compares the distribution of the J-band absolute 
magnitudes of thes e galaxies with the J-band luminosity function 
of local galaxies (Jones et al.L 120061) . A Kolmogorov-Smirnov 
test indicates a probability of 98% that this sample and the local 
luminosity function have a similar distribution. Thus our local 
sample of galaxies is representative of the My < -20.3 galaxies 
in the local Universe. Additional info rmation about this sa mple 
observed by the SDSS can be found in Fuku gita et al.l d2007l) and 
references therein. 



2. Representative samples of nearby and distant 
galaxies 

The first sample represents the local galaxies, constructed from 
SDSS survey. The second one represents the distant galaxies se- 
lected from GOODS. To gather both samples, we used a sin- 
gle selection criterion, i.e., the absolute magnitude in J band 
(My < -20.3). However for practical reasons we further re- 
stricted our sample to the galaxies possessing a good quality 
spectrum to measure [On]/13727 emission line, as well as having 



1 Intermediate MAss Galaxy Evolution Sequence is a VLT large 
program aiming at recov ering the spatially resolved kinematics of ~ 
0.4< z < 0.9 galaxies, see lYang et all 12008). 



2.2. Distant galaxy sample 

This sample is made up of the combination of two subsam- 
ples, distant star-forming galaxies and distant quiescent galax- 
ies. Those have been arbitrarily defined by their [On]/l3727 
EW, l arger or smaller than 15A, respectively (see lHammer et all 
fl997l) . In both cases, galaxies were selected with the follow- 
ing criteria: My(Afi) < -20.3 and redshift range 0.4 < z < 
0.8. All these limits were adopted to be consistent with the 
IMAGES survey. Furthermore, we selected only galaxies hav- 
ing GOODS(ACS) v2.0 images in at least 3 bands (v - 5915 A, 
i - 7697 A and z- 9103 A). 



2 This include bad quality imagery and bad pixels near the [On]/l3727 
emission line. 
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There is a relatively good correspondence between star- 
forming and quiescent galaxies defined as above with blue and 
red galaxies, respectively. Figure 3 compares rest-frame (U-B) 
colors of galaxies in th e distant sample with the limit applied by 
ICirasuolo et al.l (120071) to separate blue from red galaxies. It is, 
however, slightly imbalanced as there are more quiescent galax- 
ies in the blue cloud than star-forming ones in the red cloud. 
An optimal balance would be provided if we had chosen an 
[On]i3727 EW of about 1 1A This might slightly alter the com- 
parison with ICirasuolo et ail d2007l) in Fig. 2 (b, c, see also be- 
low). 

Distant star-forming galaxies: This sample of 52 starburst 
galaxies (EW ([OnU3727) > 15A) has been gathered by 
INeichel et all (120081) in the frame of the IMAGES survey. 
We kept only 49 galaxies, and two were rejected because 
they do not have at least 3 bands GOODS(ACS) v2.0 im- 
ages (J030238.74+000611.5 & J223256.08-603414.1) and one 
of them (J033210. 76-274234.6) ha s been recently identify as a 
quiescent galaxy dYang et al.l 120091) . Figure 2(b) shows the dis- 
tribution of the J-band absolute magnitudes of the distant star- 
forming galaxy subsample. To further verify its representative- 
ness, we over-plotted in this histogram the K-band luminosity 
function obtained by ICirasuolo et alJ (120071) for blue galaxies 
in the redshift range 0.25 < z < 0.75. The correction from 
K-band to J-band is very small in the AB system, according 
to the available values calculated in the IMAGES survey. The 
J-K values average to -0.08+0.03 and are almost independen t 
of the spectrophotometric type (see also lMannucci et al ., 2001). 
In the following we assume a J-K correction of -0.1 mag. Our 
sample follows the same distribution of the luminosity func- 
tion (Kolmogorov-Smirnov test, with a probability of 97% ), as 
illustrated in Fig. 2 (b). It confirms the Yang et al. (2008j) and 
INeichel etaD (120081) results that the IMAGES survey is repre- 
sentative of distant starbursts. 

Distant quiescent galaxies: This subsample complements 
the previous subsample, as it includes only galaxies with 
(EW( [On] ,13727) < 15A). It was selected in the GOODS 
area, using spectroscopic data collected fro m (1) observa- 
tions made by VLT/FORS2 by IMAGES (iRodrigues et al.L 
120081). (2) FORS2 observ a tion i n the frame of ESO follow- 
ups (IVanzella et all 120051 l2006l). (3) observations made by 
VLT/VIMOS dLe Fevre et aUl2005l) . For the VVDS survey, only 
secure redshift were used (with zflag= 3/4), which correspond to 
more than 95% of confidence. The first step was to detect those 
spectra with a wavelength coverage including the [On]/i3727 
line (the same criterion that is applied to the local sample galax- 
ies). We then took only those galaxies with EW([On]/i3727) < 
15A) and an absolute magnitude My < -20.3 within the redshift 
range 0.4 < z < 0.8 (same criteria of the IMAGES large pro- 
gram). Finally, we kept only the galaxies having GOODS(ACS) 
v2.0 images in at least 3 bands, leading to a sample of 94 quies- 
cent galaxies. A high fraction (77% of the sample) have spectra 
observed with VLT/FORS2 ensuring the best accuracy in mea- 
suring [On]/13727 EW. Figure 2(c) shows the corresponding dis- 
tribution of the J-band absolute magnitudes that is c ompared 
to the luminosity function from Ciras uolo et al.l d2007l) for red 
galaxies in the redshift range 0.25 < z < 0.75. The J-K mag- 
nitude correction was done in the same way as in the previous 
subsample. A Kolmogorov-Smirnov test provides a probability 
of 94%, that the quiescent galaxies sample and the correspond- 




z 



Fig. 1. Galaxies of the Fuku gita et alJ d2007l) sample in the Mj 
versus z plane. The red points are those included in the selected 
representative sub-sample of 218 galaxies with 0.02 < z < 0.03. 
This cut in redshift was determined to optimize the representa- 
tivity in terms of the LF. Curved blue line shows the limit of the 
sample that is assumed to be represented by a blue starburst with 
r P =16. 
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Fig. 3. U-B color distribution of the distant sample (blue squares: 
starbursts, red crosses: quiescent galaxies). The black line rep- 
resent s the color-magnitude relation applied by Cira suolo et al.l 
d2007l) to separate their blue and red populations (see their Eq. 
2). 

ing luminosity function in the redshift range 0.25 < z < 0.75 are 
drawn from the same distribution. Both distributions show well 
the peak that is due to massive early type galaxies. 

3. A similar morphological analysis for both local 
and distant galaxies 

3.1. Comparing the available data for the two samples 

Our main goal is to ensure that the comparison of morphologies 
is not affected by biases related to spatial samplin g, depth, and 
^-corr ections. In this respect, the detailed study of B arden et al.l 
(2008) is very encouraging. They indeed created a set of SDSS 
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Fig. 2. (a) Local galaxy sample. The black line shows the histogram of My for all 116 galaxies. Observed LF in the redshift range 
0.0025 < z < 0.2 (green line) is also plotted as a reference. This observed LF come from ! Jones et a l. (2006). (b) Distant star-forming 
galaxy sample. The black line shows the histogram of My for all 49 galaxi es. Observed LF for 0.25 < z < 0.75 (blue line) is also 
plotted as a reference. This observed LF comes from lCirasuolo et ail d2007l) for blue galaxies, (c) Distant quiescent galaxy sample. 
The black line shows the histogram of My f or all 94 galaxies. Obse rved LF for 0.25 < z < 0.75 (straight red line) is also plotted 
as a reference. This observed LF comes from lCirasuolo et ail {2007) for red galaxies. In all three plots, the vertical dashed red line 
indicates the My IMAGES limit of -20.3. 



galaxies that they have artificially redshifted to redshifts 0. 1 < 
z < 1.1, in order to compare them with GEMS and COSMOS 
gal axies. 

iBarden et al.l d2008) applied and verified whether these tem- 
plates may be compared to galaxies observed with HST/ACS. 
They did not find any significant difference in estimating ab- 
solute magnitude, half light radius, and Sersic index (see their 
Fig. 9). As they concluded, their software has no unwanted sys- 
tematic effects on measured morphological and photometrical 
parameters of simulated galaxies, and deviations from theoreti- 
cal values are well within t he uncertainties expected from their 
GALFIT dPeng et al.ll2002l) simulations. 

At z=0.65, the concordance of the sampling in physical units 
is optimal: the HST/ACS of GOODS delivers an FWHM of 
0.108 arcsec that corresponds to 0.81 kpc. At z=0.025 an SDSS 
galaxy is observed within an average FWHM=1.4 arcsec that 
corresponds to 0.74 kpc. Both pixel sizes correspond to a sam- 
pling of the PSF by 3.6 and 3.5 pixels , respe ctively. Thus the 
conditions reproduced by IBarden etaD (2008) are particularily 
appropriate for a comparative analysis of the local and distant 
samples of galaxies considered in this study. 

We thus have to consider two additional effects, i.e., the com- 
parison of the optical depth of the two databases and the k- 
correction effects. Concerning the latter, there is excellent agree- 
ment between the rest-frame band of distant galaxies to those 
of local galaxies, assuming z=0.65 for the median value of dis- 
tant galaxies (see Table 1). Table 2 summarizes the observational 
conditions resulting from both SDSS and GOODS. We are in- 
terested in verifying whether the analysis of distant galaxies is 
based on data as deep as or deeper than those used for local 
galaxies: 
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It turns out that, after accounting for cosmological dimming 
and ^-correction (the last term in Eq. 1), the GOODS imaging 
is deeper than that of SDSS, by 0.52, 0.08, and 1.02 magnitude, 



Survey 




u band 


g band 


r band 


i band 


z band 


SDSS 




3551 A 


4686 A 


6165 A 


7481 A 


8931 A 




B band 


Vband 


i band 


z band 






GOODS 


4312 A 


5915 A 


7697 A 


9103 A 






rest-frame 


2582 A 


3542 A 


4609 A 


5451 A 







Table 1. Wavelength comparison of different bands in SDSS and 
GOODS(ACS). 



in rest-frame u, g, and r bands, respectively. Possible effects re- 
lated to depth will be discussed later in the discussi on section- 
Beside s this effect, we confirm the conclusion of B arden et al.l 
(120081) that there should not be systematic effects on measuring 
morphological and photometrical parameters from one sample 
to the next. 



3.2. Light profile analysis 

Half-light radius and bulge-to-total ratio (B/T) are derived from 
the surface brightness profile analysis. In this work, instead of 
fitting simultaneously all neighboring objects within 5.0 arcsec 
of the target, they were masked to avoid any flux contamina- 
tion. To measure the half-light radius, we have developed our 
own IDL procedure, which allows us to visually analyze differ- 
ent profiles of the galaxy at the same time: (1) the flux profile, 
(2) the magnitude profile, and (3) the flux content within ellipses 
of different radius with a step of one pixel. The first two let us 
determine a reasonable sky value and the last one gives us the 
half light-radius when a real "plateau" is visualized. 

Each galaxy is decomposed as a combination of a bulge and 
a disk. This decomposition is a two-dimensional modeling based 
on light (flux) distributions. We then modeled each galaxy light 
profile with a combination of two Sersic laws. The first one 
would represent the bulge with a Sersic index that is left free, 
and the second one would represent the disk with a Sersic index 
equal to one. B/T is derived from the GALFIT decomposition of 
the galaxy. 

Half-light radius measurements and GALFIT simulations 
were done at rest frame r band, i.e., at z band for galaxies in 
the distant sample. 
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SDSS 


GOODS ACS 


D=telescope diameter (m) 


2.5 


2.4 


Band 


u 


g 


r 


B 


V 


i 


z 


T=Expo-time (s) 


53.907456 


53.907456 


53.907456 


7200.00 


5450.00 


7028.00 


18232.00 


B=sky background (mag) 


22.15 


21.85 


20.85 


23.43 


22.74 


22.72 


22.36 


Filter FWHM (A) 


567.00 


1387.00 


1373.00 


728.95 


1565.50 


1017.40 


1269.10 


Filter range (A) 


-1000.00 


-1800.00 


-1500.00 


8780.00 


2570.00 


1910.00 


>3080.00 



Table 2. Observational parameters for both SDSS and GOODS(ACS) imaging. 



3.3. Color images and maps 

Color is an important information in our morphological classifi- 
cation process. We have two samples of galaxies (local and dis- 
tant samples) and we want to compare them. It is needed to use 
same rest-frame colors in both samples to construct color maps 
(see Table 1). Therefore, in the case of two-color maps, we use 
u-r bands for the local sample and v-z bands for the distant sam- 
ple. We do the same for three-bands color images (u-g-r and v-i-z 
bands, respectively). 

Three-band color images are useful for the morphological 
classification for examining the small-scale structures of the 
galaxy. However, we need color maps in which color of individ- 
ual features can be measured. To do so, we substract, pixel by 
pixel, the magnitude in two observed bands using an algorithm 
that all ows estimation of c olors and their uncertainties (see de- 
tails in Zheng et all 120051) . A color map of any substructure in 
the galaxy can then be compared with models of stellar popula- 
tion synthesis. They can thus be useful for estimating the color 
of a bar or identifying dusty or star-forming regions. 

3.4. Morphological classification 

To make the morphological classification a reproducible process 
and to reduce its inherent subjectivity, we constructed a pseudo- 
auto matic decision tree (d escribed in Fig. 4). It was adapted 
from iNeichel et alj (120081) to account for SO and E galaxies. 
The decision tree takes all the available quantities into account 
such as B/T, half light radius, GALFIT parameters, GALFIT 
model and residual images, error images, disk-bulge-galaxy 
profile, two color-maps, and three-color images. To reduce 
the residual subjectivity, the morphological classification of 
each galaxy of the samples was done by three of us (RD, FH, 
YY) following the decision tree. Then, our classification was 
compared. Because of the relative simplicity of the decision 
tree, the agreement between the various classifications was 
excellent for both samples of local and distant galaxies. 

A few changes to the INeichel et al.l (120081) methodology 
were made in the present morphological classification. First, 
we now formally account for the agreement between bulge and 
disk centers. Indeed, we have noticed that, in some cases, the 
fit bulge+disk of galaxies seems to be good, but in the three- 
color image or in the residual fit image we see disagreement 
between bulge and disk centers. Galaxies with a disagreement 
in their centroid larger than 3 pixels are considered as peculiar 
galaxies ( see decision tree). Se cond, and conversely to what was 
chosen in iNeicheTe t al. (2008), we classify as compact all ob- 
jects with a half-light radius below 1 kpc. In previous studies of 
z ~0.6 galaxies, the limit to defining compact galaxies has been 
taken to be 3 kpc dMelbourne etall [2005; Zh eng et all [2005; 
lHammer et al.L 1200 lb because in all those works the spatial res- 
olution of WFPC2 images was significantly higher than the 



0".108 FWHM provided by ACS/HST images for most (90%) 
galaxies in the distant sample. We then have access to a bet- 
ter spatial resolution. Thus, the compactness limit is 1 kpc, and 
only such and extreme compactness is limiting our morphologi- 
cal classification. 

To summarize, we distinguish between the following classes: 

1. E: elliptical galaxies, with a B/T ratio between 0.8 and 1.0. 

2. SO: lenticular galaxies. There is no presence of regular struc- 
tures (arms) and the bulge is redder than the disk. A highly 
symmetric disk, as well as the match of bulge and disk cen- 
ter, is observed. The B/T ratio has a value between 0.5 and 
0.8. 

3. Sp: spiral disk galaxies. This class is characterized by a bulge 
redder than the disk. Regular arms and a highly symmetric 
disk are also present. Bulge and disk centers coincide, and 
the B/T ratio is smaller than 0.5. 

4. Pec: peculiar galaxies. The main characteristic of this class 
is the presence of asymmetrical features. However, we may 
divide this class into four subclasses: possible mergers 
(Pec/M), which are objects with irregularities that could be 
associated to merger/interaction events; tadpole-like (Pec/T), 
which are objects showing a knot at one end plus an extended 
tail; irregulars (Pec/Irr) for objects with asymmetric irregu- 
larities that cannot be associated with arms, objects with a 
center (bulge) bluer than the rest, or objects with a disagree- 
ment between bulge and disk centers larger than 3 pixels; and 
compact (Pec/C), which are galaxies with a half-light radius 
smaller than 1.0 kpc. 

According to this decision tree, we have classified the 
259 galaxies that belong to our two samples. Galaxies from 
INeichel et al.l d2008l) have kept their morphological classifica- 
tion, except for those classified as compact by them but that are 
not in the present study. 



4. Morphological classification: results 

Table 3 summarizes the fraction of E, SO, spirals, and peculiar 
galaxies. Peculiar galaxies are also subdivided into 4 subcate- 
gories as described in Sect. 3.4. 

4.1. Morphological classification of local sample 

By applying our morphological classification method to the lo- 
cal sample, we confirm that the local Universe is dominated 
by spirals galaxies that represent 72% of the population, and 
that the fraction of E/S0 is 18%, while 10% of galaxies have 
peculiar morphologies. This confirms the fraction of spirals in 
local galaxies as found by Ham mer et al.l (120051) on the ba- 
sis of the morp hological determined luminosity function of 
iNakamura et al.l d2004l) . The number of compact galaxies in the 
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Fig. 4. Semi-automatic decision tree used in the morphological classification process. Each step considers a simple and unique 
criterion. We chose a very conservative method, which takes into account the well known morphologies of local galaxies that 
populate the Hubble sequence. 



local sample is very small (2%), which proves that spatial res- 
olution can not significantly alter our classification. Among the 
116 galaxies of the local sample we also find 96 (83%) quies- 
cent galaxies (EW([On]^3727]) < 15 A) and 20 (17%) st arburst 
galaxies (EW([On]/l3727) > 15A). lHammer et alj dl997l) . using 
a completely different method, find the same percentage of star- 
burst galaxies which have same definition than in the present 
paper (see their Fig. 3 and their Sect. 3.1). 



5. Building up the past Hubble sequence 

5. 1 . Morphological evolution by comparing distant and local 
galaxies 

The result s shown in Ta b le 3 a re in a very good agreement with 
those of van den Berghl (l2002l) who found that the fraction of 
peculiar galaxies increases from 12% to 46% from z=0 to z=0.6- 
0.8. A simple examination of Table 3 reveals that during the past 
6 Gyrs: 



4.2. Morphological classification of distant sample 



As assumed bv lHammer et all (120051) or bv lNeichel et al.l (120081) . 
one would have expected a small fraction of peculiar morpholo- 
gies in the sample of quiescent galaxies. In contrast to that ex- 
pectation, we find a quite surprisingly high fraction of peculiar 
galaxies, 25%. Peculiar galaxies are also the more common type 
(69%) in the starburst distant sample, in which no elliptical or 
lenticular galaxies are found. Gathering the whole sample of 
distant galaxies, we find that 52% of them have peculiar mor- 
phologies, while only 31% have spiral morphologies, with the 
rest either E or SO. Gathering distant starbursts and quiescent 
galaxies in a s i ngle r epresentative sample was done following 
lHammer et al.l (1 19971) : galaxies with significant line emission 
(EW([On]^3727) > 15A) represent 60% of intermediate red- 
shift (0.4 < z < 1.0) galaxies, while quiescent galaxies represent 
40%. Such a distribution is fa irly well confirm ed by the fraction 
of red and blue galaxies from lCirasuolo et all d2007l) in the same 
redshift range. 



- the fraction of early type galaxies (E/S0) has not evolved, 
suggesting this population was mostly in place at that 
epoch, consis t ently with conclusions o f other studies 
dJimenez et all 120071: iBernardi et alll2006h : 

- the fraction of regular spiral galaxies has increased by a 
factor 2.3 and this result supersedes former claims (e.g., 
iLillv et all . 1 19981) that were based on much lower data qual- 
ity; 

- the fraction of peculiar galaxies has decreased by a factor 
5, and this represents a very strong evolution, as half of the 
galaxy population consisted of peculiar galaxies 6 Gyrs ago. 

At first glance, the number evolution of spirals is balanced 
by the number evolution of peculiar galaxies. Be cause peculiar 
galax ies generally show anomalous kinematics dNeichel et all 
2008), this suggests that peculiar galaxies have to be transformed 
in one way or another into rotating spirals that we observe in 
large numbers today. Such a transformation has to be under- 
stood and reported by the different models of galaxy evolution. 
However, before taking that step forward, we need to verify 
whether this representation of the past Hubble sequence can be 
causally linked to present-day galaxies. 
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Type 


Total (%) 


Quiescent (%) 


Starburst (%) 


Total(%) 


Quiescent (%) 


Starburst (%) 


E 


3±1 


3±2 


0+0 


4±1 


11±3 


0+0 


SO 


1 ^-i-A 
1 J±4 


1 /I -uA 
1"4±H- 


ZU± 1U 




JJ±0 


u±u 


Spiral 


72+8 


76±10 


55+17 


31+7 


31+6 


31±8 


Peculiar 11 : 


10±3 


7±3 


25+11 


52±9 


25±5 


69±12 


P/Irr 


4±2 


2+1 


15±9 


26±7 


21+5 


29±8 


P/Tad 


0±0 


0+0 


0±0 


6±3 


0±0 


10±5 


P/Mer 


4±2 


4±2 


5±5 


20±6 


4±2 


30+8 


P/C 


2±1 


1±1 


5±5 


0+0 


0+0 


0±0 



Table 3. Fraction of the different morphological types for local and distant samples. a Fractions in the peculiar class are simply 
addition of the fraction of peculiar subclasses. 



5.2. Have we established a past Hubble sequence? 

In this paper we have tried to keep a homogeneous methodol- 
ogy of classifying both distant and local galaxies on the basis 
of a similar observational set-up. The only difference that could 
bias our result is the greater optical depth -by 0.08 to 1 .0 mag- 
reached for distant galaxies when compared to local galaxies. To 
verify whether this might affect our results in Table 3, we have 
reproduced similar conditions at z=0.65 than at z-0 by degrad- 
ing the signal to noise of GOODS observations. It results that 
small variations of optical depth cannot affect our results in Table 
3: all the galaxies which have been classified as peculiar remain 
in that category. This is not surprising because SDSS imaging 
is deep enough to p rovide a fairly robust classification of galax y 
morphology (e.g jNakamura et all 120041; IFukugita et al.Ll2007h . 
Figure 5 shows the two Hubble sequences that have been con- 
structed using both local and distant samples. 

The main drawbacks in establishing such a sequence are, 
however, related to our fundamental ignorance of the various 
effects on galaxy morphology during the past 6 Gyrs. As dis- 
cussed in Sect. 1, the possible effects on either the number den- 
sity of galaxies or the validity of a common absolute magni- 
tude limit include galaxy merging and stellar population evo- 
lution. Their possible impacts are examined in Appendix A. It 
turns out that the past Hubble sequence displayed in Fig. 5 pro- 
vides a reasonable snapshot of the past history of galaxies pop- 
ulating the Hubble sequence today. The link between the two 
Hubble sequences is marginally affected by very recent merg- 
ers or by stellar population evolution. Some very specific star 
formation histories could, however, complicate the link: for ex- 
ample a galaxy forming all its stars at z=l in a so-called mono- 
lithic collapse would show an evolution of its M s/e /; flr /L^ ra- 
tio by about 0.4 dex according to population synthesis model. 
Such objects, if they exist, would be found in slightly greater 
numbers in the past Hubble sequence than in the present-day 
one. Besides this, the proposed link between past and present 
sequen ce is also well adapted to various models of galaxy evo- 
lution (Bruzua f& Charloll2003h . For example, let us consider a 
secular evolution model in which a galaxy gradu ally accretes gas 
from the intergalactic medium. Observations (Ha mmer et al.L 
12005b iBeil et all [2005 ) suggest a stellar mass increase by 0.3 
dex since z=l, or 0.15-0.2 dex since z=0.65. By selecting local 
and distant galaxies through Mj(AB)< -20.3, one would catch 
both progenitor and descendant because the observed evolution 
of the M s teiiar I ' Lk ratio is 0.15 dex. 

5.3. A few results and possible applications 

There are a few results displayed in Table 3 and in Fig. 5 that 
deserve some comment, as they appear to be unexpected: 



1. The fraction of 10% of peculiar galaxies in the local Hubble 
sequence is not negligible, given that these galaxies have 
stellar masses well above 10 10 M . Indeed. IFukugita et alj 
(120071) find only 1 % of galaxies with peculiar morphologies. 
This may be related to the not using the color information 
and thus not accounting for some galaxies (at a level of ~ 3% 
in our sample) having a blue "bulge" surrounded by a redder 
disk. A part of this discrepancy (an additional 2-3%) comes 
from our additional criterion that the bulge and disk centroids 
should coincide to classify a galaxy as a spiral. Otherwise, 
local peculiar galaxies seem to be dominated by mergers or 
compact galaxies. 

2. We also find quite a high fraction ( 25% ) of peculiar mor- 
phologies in the sample of quiescent and distant galaxies. 
This is not mainly caused by the effects des cribed above. 
This r esult slightly affect t he conclusio ns of iNeichel et al.l 
(2008) and HammeretalJ (|2005j [2009), who assumed that 
besides E/S0, all quiescent galaxies have to be spirals. This 
simply reduces further the fraction of regular, rotating disks 
in the past Hubble diagram whose fraction is 2.3 times lower 
than in the present-day galaxies. 

Figure 5 presents an accurate representation of Hubble se- 
quences at two different epochs for which galaxy morphologies 
have been studied in an homogeneous way and with a minimum 
of selection effects. This can also be used to test some scaling re- 
lations betwe en fundamental par ame ters. For example, this ha s 
been done by iFlores et"al . (2006) and lPuech et all (120081 120091) . 
who show that the main evolution of the Tully-Ficher relation 
is resulting from a much larger scatter from peculiar galaxies. 
Figure 2 demonstrates that this strong evolution discovered by 
IMAGES is unaffected by selection biases. 

Similarly, one can test the scaling rela tions of early type 
galaxies. For examp le, it has been shown dPahre et all [l998; 
iBernardi et alll2003l and references therein) that bulge luminos- 
ity correlates with the bulge effective radius, which is indeed one 
projection of the fundamental plane of elliptical galaxies. Figure 
6 presents the absolute r-band magnitude of the bulge against the 
bulge radius for the two subsamples of E/S0 galaxies. Although 
the number of objects is quite small, there is a correlation be- 
tween the luminosity and bulge radius. Using the past Hubble 
sequence, we may verify whether it evolves or not. Figure 6 (bot- 
tom panel) shows that no evolution is measured, when account- 
ing for passive evolution of the luminosity (about 0.59 mag), for 
stellar populations assumed to be formed ~ 12 Gyr ago. 

6. Conclusion 

We have established a first approximation of what would be 
the progenitors, 6 Gyrs ago, of the galaxies of the present-day 
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Fig. 5. Present-day Hubble sequence derived from the local sample and past Hubble sequence derived from the distant sample. Each 
stamp represents approximately 5% of the galaxy population. 
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Fig. 6. r absolute magnitude of the bulge vs bulge radius of E 
(red crosses) and of SO (yellow squares) galaxies. Top: Local 
sample for which the soli d line represent s the b est fit, while the 
dotted line comes from the Bernardi et al . (2003) study of ellipti- 
cal galaxies. Middle: Same for distant E and SO galaxies with the 
dashed line showing the best fit of the relation. There is approx- 
iamtely 0.5 magnitude between the local (solid line) and distant 
best-fit of the M r -Rb„ig e relation. Bottom: Same as above but af- 
ter correcting for a passive evolution of the stellar population. 
We assume a single star formation burst, 12 Gyrs ago. Other star 
formation histories could be considered, but this indicates that 
the M r -Rbui ge is unlikely to be a strongly evolving relationship. 



Hubble sequence. Using a simple and single criterion for select- 
ing both distant and local galaxies, we showed that one may link 
them together quite robustly. There are certainly uncertainties 
caused by further, very recent merging events, as well as by var- 
ious star formation histories, but they are not likely to affect more 
than 10% of the galaxies that may be linked together. This is 
lower than the error due to Poisson number statistics in samples 
that slightly pass 100 galaxies. Our method classifying galaxy 
morphologies is similar to that adopted for local galaxies. It uses 
a semi-automatic decision tree that is applied in a consistent way 
for both distant and local galaxies. The errors due to misclassifi- 
cation are thus unlikely to be more than a few percent. This leads 
to: 

1 . E/S0 galaxy populations show no evidence of number evolu- 
tion during the past 6 Gyrs; 

2. Slightly more than half of the distant galaxies have peculiar 
morphologies, whic h is likely associated to anomalous kine- 
matics according to lNeichel et alJ (120081) : 

3. The fraction of regular spiral was 2.3 times lower 6 Gyrs ago 
than in the present epoch; 



4. Almost all the evolution is caused by the transformation 
of galaxies with peculiar morphologies into regular spiral 
galaxies at present epoch. 

The past Hubble sequence that is established in this study 
could be a useful tool for further study of galaxy evolution or 
transformation. It can be used to derive the evolution of fun- 
damental planes for both spirals and bulges. It is, however, un- 
clear whether this method can be extended to longer look-back 
times, at epochs during which mergers could be a predominant 
mechanism to form galaxies and their stellar populations. There 
is, however, a need to improve the statistics, since this paper is 
based on numbers slightly over 100 for both distant and local 
samples, i.e., fewer galaxies in the different morphological types. 
Cosmic variance effects also have to be investi gated since they 
could easily affect the GOODS South field (see Ravikumar et all 
2007, and references therein). Further selection of SDSS galax- 
ies and galaxies from GOODS South and North or from deep 
galaxy surveys with similar depth are feasible. More than half 
of the galaxies in the past Hubble sequence have peculiar mor- 
phologies. Their transformation to regular spirals in the present- 
day Hubble sequence is a big challenge that should be addressed 
by current scenarios of galaxy evolution and formation. 
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Appendix A: Possible biases and selection effects 
in linking the two Hubble sequences 

A. 1. Effects related to merging 

Merging naturally affects the evolutionary link between distant 
and present-day galaxies, simply by a radical suppression of one 
of the merging galaxies (see bottom of Fig. A.l). Another and 
more subtle effect would be caused by two galaxies below our 
limit (e.g. with Mj(AB) > -20.3) merging and their by-product 
reaching our limit (Mj(AB) <-20.3) at z=0 (see top of Fig. A.l). 



z= 0.65 


Z= 




Mj > -20.3 + 1 galaxy 

s*s~s 

1 Mj < -20.3 
Mj > -20.3 


Chance < 4% 


m j 5 - 20 - 3 - 1 galaxy 

1 Mj < -20.3 
Mj < -20.3 


Chance < 2% 



Fig. A.l. Two possible merging effects that could affect the link 
between distant and present-day galaxies. Bottom: two galax- 
ies with Mj(AB) < -20.3 merge after z=0.65 suppressing one 
galaxy with Mj(AB) <-20.3 at z=0.0; Top: two galaxies with 
Mj(AB) > -20.3 merge after z=0.65 giving one galaxy with 
Mj(AB) <-20.3 at z=0.0. 



To evaluate the fraction of z=0.65 galaxies that will merge 
during the last 6 Gyrs, we may consider the two studies 
providing the largest and the smallest evolutio n of the merger 
rates. Assuming an evolution rate as (1 + z) 2 ' 1 dLe Fevre et all 
2000), one can derive that 29% of z=0.65 have to experience 
a maj or merger during the last 6 Gyrs. Similarly, lLotz et all 
(2006) estimate that between 33% and 66% of L B > 0.4x L B 
galaxies had a major merger since z=l.l, with an evolution rate 
of (1 + z) 1 ' 12 . Combining the above estimates implies that 15 



to 30% of z=0.65 galaxies could experience a major merger 
later on. This number is not very accurate and depends on 
assumptions about the selection filter, as well as on the exp ected 
time t for a galaxy pair to merge. iRawat et"aTI (|2008) have 
estimated the merger rate using a quite similar method than ours 
to pre-select galaxy pairs, i.e., within Mj(AB) < -19.5 galaxies 
selected in the GOODS field. They do find an evolution follow- 
ing 0.014x(l + z) 2 1 ■ For t=0.35 and 0.5 Gyr, one calculates that 
21 and 15% is the fraction of z=0.65 galaxies that are expected 
to experiment a major merger before reaching z=0, respectively. 

Major mergers may dramatically affect the appearance and 
dynamical properties of galaxies. Here we are interested in link- 
ing distant progenitors with their local descendants. In the fol- 
lowing we assume that the actual limit in linking past to present 
Hubble sequence is realized when a galaxy is sim ply suppressed 
by be ing accreted by a more massive galaxy. [Hammer et alj 
(2009), assuming that all objects could be modeled as a merger 
or a merger remnant, have derived the corresponding merger 
mass ratio (see their Fig. 6, bottom). The probability for an 
Mj(AB) < -20.3 to encounter a more massive galaxy is about 
7%. Such a value proves the scarcing of events involving a 
galaxy more massive than the observed one, since those are rarer 
due to the exponential drop of the mass function towards the 
massive end. Thus the chance for a galaxy to be suppressed in 
the local Hubble sequence by a very recent merger, would be 
0.21x0.07, hence less than 2%. 

We also have to account for the possibility that two Mj(AB) 
> -20.3 galaxies could merge during the past 6 Gyrs and reach 
the local Hubble sequence of Mj(AB) < -20.3 galaxies (see top 
of Fig. A.l.). The probability of such events is similarly small. 
To illustrate this, let us assume two progenitors having the same 
J luminosity: only merging galaxies with Mj(AB) < -19.6 can 
be progenitors of an Mj(AB) < -20.3 galaxy at z=0. From 
Fig. 2 such galaxies represent less than 20% of the Mj(AB) 
< -20.3 sample. By integrating all possible major mergers 
in a similar way, it leaves less than 4% of z=0.65 galaxies 
that may enter the sample due to a recent major merger. We 
conclude that very recent mergers are unlikely to affect the link 
between past and present Hubble sequences at a level over ~ 5%. 



A.2. Effects related to stellar population evolution 

Our selection is based on the absolute J-band magnitude that 
is often considered as a proxy for stellar mass. There is also 
very good correspondence between absolute K-band and J-band 
magnitudes (see Sect. 2.3). The stellar mass to K-band light 
ratio is evolving quite significantly with the redshift, which is 
simply related to the evolut ion in co l or or in star formation. 
For example, ap plying the | Bell et al] (l2003h correction to the 
M ste u ar /Lk ratio, iPuech et aD d2008l) find an average evolution 
of 0.15 dex when comparing local spiral galaxies to distant 
emission line galaxies. In other words, by selecting more and 
more starbursts in the past Universe, we pre-select galaxies 
with lower stellar masses than those in the local Universe. 
Such an effect might be non negligible because there are 9% 
of distant starbursts that have -20.675< Mj(AB) < -20.3. 
Interestingly a part of this bias is compensated for when 
considerin g the baryonic mas s rather than the stellar mass (see 
Fig. A. 2). [Hammer et alj (120091) find that distant starbursts have 
a gas fraction averaging 31%. This is based on assuming the 
Kennicutt law that links surface density of star formation rate to 
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that of gas mass. Considering baryonic masfl we do find that 
the larger amount of gas in distant starbursts add a shift of -0.16 
dex on to Mbaryonic/LK ratio, thus compensating quite well for 
the stellar population evolution. Of course this does not mean 
that, for any evolutionary s cheme, such a compe nsation may 
be realized as, for example, iMaraston et al.l d2006l) claim for a 
longer evolution of the M ste \\ ar jLK ratio. 



3 We asssume the gas fraction to be 10% for local galaxies. 
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Fig.A.2. Stellar and baryonic mass of a galaxy . Left: galaxy with Mj(AB) < -20.3 in the local Universe; Right: galaxy with 
Mj(AB) < -20.3 at z=0.65. 



